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I. INTRODUCTION 

One of the most exciting developments in cold atom 
experiments is the ability to emulate the Hamiltonians 
of charged particles in magnetic fields [T and electrons 
with spin-orbit coupling These techniques allow one 
to tune the dispersion in complicated spatially dependent 
ways. Here we show the that resulting single particle 
quantum mechanics is profoundly different than what we 
are used to. For example, theorems about the number of 
nodes in the ground state [5] no longer apply, and by tun- 
ing experimental parameters one can change the ground 
state in a double well from symmetric to antisymmetric, 
and back again. 

In the experiments of Lin et al. [2 HI IH [S], Rubid- 
ium atoms in the F=l hyperfine manifold interact with 
two co-propagating lasers. These lasers drive Raman 
transitions between the three magnetic hyperfine states 
771 = —1,0, 1. Including the effect of the quadratic Zee- 



man field, the Hamiltonian in the rotating frame is 

Hs = —k +-)l + 5S, + hwg\0 
+ ^Sa; cos {2kLx) ^Sy sin {2kLx) 

The matrices I, Sx,y^z sue the 3x3 identity and spin ma- 
trices in the basis (m = 1, 0, —1), h is Planck's constant 
divided by 2tt, m is the atomic mass, 6 is the effective 
detuning between states m = and m = —1, d + huiq is 
the effective detuning between states m = 1 and m = 0, 
^r/K is the Rabi frequency of the Raman lasers, and Ul 
is the recoil from the Raman lasers. In the experiment 
LJq is tuned via the quadratic Zeeman effect. 

In the limit where hujq is large, the m — 1 state is far 
off resonance and decouples. The resulting energies are 



(2) 

We will solely be concerned with the lowest energy band E_{k). Within a semiclassical treatment, where an external 
potential varies slowly, the low energy Hamiltonian is formally 



2m dx^ 



2m ox 



(3) 



where H, = flii/\/2. Throughout the rest of this work wc 
focus on the case where 6 = 0. 

In sections |III| through |Vj we study the single particle 
eigenstates of Eq. ^ for a range of archetypical poten- 
tials, and compare their properties with those of standard 
quantum mechanics. The results are unexpectedly rich. 
We explain how to reveal the nonconventional features in 
experiments. Section |VI| addresses the limits of validity 
of Eq. ^. 



There are two key properties of the dispersion in 
Eq. ([3|. First, the dispersion has two degenerate min- 
ima. Typically this results in a ground state wavefunc- 
tion which oscillates in space. Second, the dispersion is 
anharmonic. Some of the theorems in quantum mechan- 
ics (such as the non-existence of nodes in the ground 
state) are based on having a purely harmonic dispersion. 
In appropriately tailored potentials, an anharmonic dis- 
persion with a single minimum can even lead to ground- 
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state nodes. To draw out the role of these features, we 
consider the foUowing potentials: infinite square well, 
harmonic oscillator, double well. Anharmonic dispersion 
are ubiquitous in lattice systems, but we are unaware of 
an analogous study. Lattice systems have the additional 
feature that momentum space is periodic. 



II. TECHNIQUES 

We adimensionalize Eq. ([3 1, scaling all energies by 
El = h?k\/2m and lengths by fc^^. We restrict our- 
self to one dimension (ID), assuming that a tight trap 
has frozen out motion perpendicular to the x direction. 
Extending the discussion to the three-dimensional case 
is straightforward, but the interesting results already ap- 
pear in ID. 

In terms of the dimensionless variable y — k^x, the 
dimensionless Hamiltonian becomes 



n = 




+ V{y) (4) 



where n = VL/El, and V (y) = V{y/kL)/EL. 

To numerically study Eq. Q, we discretize space, and 
write the operator as a matrix, using finite difference 
approximations of various orders. We then numerically 
calculate the square root by standard algorithms. We 
verify that our results are independent of the discretiza- 
tion grid and the order of our approximation. 



III. INFINITE SQUARE WELL 



wavevector fit into the well. When an odd number fits 
best, the antisymmetric state has lower energy, other- 
wise the symmetric state wins. In Fig. [2] we show the 
two non-zero wavevectors of the ground-state inside the 
well. There are exactly two wavevectors as the equation 
E- (fc) — E can be manipulated to make a quadratic 
equation. This quadratic will have two real roots when 
17 < 4 and -1 - ^2/16 < E/El < -n/2. 

In Fig. [3] we give the number of nodes in the ground 
state as a function of the well width. The contrast with 
usual quantum mechanics, where the ground state has no 
nodes, is dramatic. 




FIG. 1: (color online) Energies of the lowest symmetric (solid 
blue line) and antisymmetric (dashed red line) eigenstates of 
Eq. (|3| taking Q — 2El, 3 = and treating V as an infinite 
square well of width L. Comparing to Fig. [3] one can see 
the states interchange as the ground state when it becomes 
energetically advantageous to add another node. 



The simplest potential to investigate is the infinite 
square well. This is most easily defined by taking the 
limit Vq — >■ oo of the finite square well 

It should be emphasized that even in the limit Vq — > oo, 
the boundary condition at the edge of the well is not sim- 
ply that the wavefunction vanishes, otherwise the opera- 
tor in Eq. Q is not self-adjoint. If one discretizes space 
as described in Sec. [ill and maintains a finite but large Vq, 
one automatically produces a self-adjoint Hamiltonian. 

The eigenstates can be classified by the number of 
nodes they possess. In Fig. [l] we show the energies of 
the lowest two eigenstates as a function of the width of 
the potential well. As one increases the well width, the 
lowest symmetric and antisymmetric state take turns be- 
ing ground state. These crossings can be understood 
by noting that there is a preferred wave-vector in the 
problem. In free space the lowest energy state has wave- 
vector fcinin = kL\Jl — {ft/4,EL)'^. As One increases the 
size of the well, different numbers of half-waves of this 
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FIG. 2: (color online) Momenta in the ground state of the 
square well as a function of well size (see Fig. [l] for param- 
eters). The solid lines show the Fourier components of the 
ground state while the dashed line between them is kmin/kL- 
The dotted vertical lines denote positions where the symme- 
try of the ground state changes. 
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FIG. 3: The number of nodes in the ground state of the square 
well as a function of well size (See Fig. [l] for parameters). 



IV. HARMONIC POTENTIAL 



of two peaks, centered at fcmin and — fcmin- The width of 
these peaks scales as 1/yo- Thus tp{y) = A{y) cos(fcminy), 
where A{y) is a smooth function that falls off on a length 
scale of order yo . 
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There are two reasons to study the harmonic poten- 
tial. First, the harmonic oscillator is one of the paradig- 
matic examples of quantum mechanics. Second, we will 
be able to get further insight into the structure of the 
ground state by considering a canonical transformation 
which switches position and momentum. Our main re- 
sults, illustrated in Figs. |4] and [5] are qualitatively similar 
to those of the infinite square well. The main differences 
are: (1) The harmonic potential favors states which have 
a higher amplitude at the center, resulting in a symmetric 
ground state. (2) One cannot readily define the "number 
of nodes" , as the wavefunction is spread over an infinite 
domain, and generically oscillates an infinite number of 
times. One can, however, consider analogous measures, 
such as how many nodes lie within a fixed number of 
oscillator lengths. 

We consider a potential of the form 



Viy) 



{yo/2r 



(6) 



and numerically calculate the eigenstates as before. 

Figure |4] shows the energy of the two lowest energy 
states as a function of yo- Unlike the square well, there 
are no level crossings, instead the two states simply 
asymptotically approach one-another. Figure [5] shows 
a density density plot of the ground state wavefunction 
as a function of y for a range of j/o- As one sees, the 
locations of the nodes are set by the characteristic wave- 
vector fcminj while a broader envelope is determined by 
the width of the potential yg. 

The simplest way to understand these results is to note 
that a Canonical transformation y — > —idx, —idy — > x 
converts this to a standard quantum mechanics prob- 
lem with a quadratic dispersion and a double well po- 
tential. Thus the the Fourier transform of the wavefunc- 
tion, tjj{k) = J dy e~^^^'ip{y), is the real-space wavefunc- 
tion of a conventional double well. That is, consists 



FIG. 4: (color online) Energies of the lowest symmetric 
(solid blue line) and antisymmetric (dashed red line) eigen- 
states of Eq. S taking Q = 2El, 5 = 0, and taking 
V = 2"^ EL{kLx) /j/o, corresponding to a harmonic potential 
with characteristic length xq = yo/kL- 
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FIG. 5: (Color online) A density plot of the wavefunction of 
the ground state of the simple harmonic oscillator as a func- 
tion of its width. Red and blue denote positive and negative 
values of ip, and darker colors correspond to larger values. 
Parameters are given in Fig. [4] 
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V. DOUBLE WELL 

Our final potential is the double well, 

-I 2 



(7) 



This is particularly interesting, as now we may have both 
a double well in momentum space and a double well in 
real space. Either of these double-wells can take the 
ground state from symmetric to antisymmetric. Figure [6] 
shows a cut through parameter space that illustrates this 
principle. 

Considering only changes to the real-space potential, 
there are several different scenarios which lead to an an- 
tisymmetric ground state. First, as in Sec. |III[ changing 
the width of the real-space well changes the number of 
half-wavelengths of wavevector /cmin that fit. When an 
even number of half-wavelengths is optimal, the ground 
state is antisymmetric. Equivalently, changing /cmin while 
fixing the real-space potential, can drive a transition. 
Second, the bump in the double-well potential, favors 
wavefunctions which have a node in the center. In regular 
quantum mechanics, this effect never drives the energy 
of the antisymmetric state below that of the symmetric 
state. Here, with the non-quadratic dispersion, one can 
however find a level crossing. By transforming k ^ x 
and X ^ k one can repeat these arguments in Fourier 
space. 




FIG. 6: The symmetry of the ground state of the double 
double well, Eq. ([7|. Increasing the ordinate, /3, increases the 
size of the barrier in the real space potential, while increasing 
the abscissa, Q/El, decreases the size of the barrier in the 
double-well momentum space dispersion. Light areas have a 
symmetric ground state and dark areas antisymmetric. Here 
we take yo — 9, S = 0. 



VI. VALIDITY OF THE SEMICLASSICAL 
APPROXIMATION 

Our discussion so far centered on the dispersion curve 
Eq. ([2]) and the resulting Hamiltonian Eq. ([3| . This form 
is achieved by applying a unitary transformation S to the 
Hamiltonian found in Eq. ([T]) , after neglecting the decou- 
pled off-resonance state. However, this transformation is 
a function of momentum, S — S (k) and so it does not 
commute with the physical potential. Under this trans- 
formation, the operator in Eq. Q representing the po- 
tential is 



'''^phys 



(8) 



where Vphys = Vphys(a;) is the physical trap potential. In 
our discussion we have neglected the second term on the 
right hand side of Eq. Q . Here we consider the limits of 
validity of this approximation. 

As is clear from dimensional analysis, the corrections to 
our approximation will involve terms such as k2^V'{x). 
If the characteristic scale of the changes in V^hys are 
large compared to l/k^, then these corrections can be 
neglected. In our dimensionless units, this requires the 
potential V (y) to only change on a length-scale large 
compared to unity. In all sections, we investigated po- 
tentials of this form, and hence we expect our results 



to be robust. Numerical investigation of the full spinor 
Hamiltonian, as shown in Fig. [7j confirms that Eq. Js]) 
quantitatively captures the low energy physics of Eq. (IsF. 
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FIG. 7: (color online) The ratio of the energies E and 
iSphya calculated respectively from Eq. ([3| and Eq [l] with 
V = 2'^EL{kLxf/y^, n = 2ELU)q = 00, 5 = 0. Blue solid Hue: 
ground states; Red dashed line: first excited states; Black 
dotted line: excitation energy. 
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VII. CONCLUSIONS AND OUTLOOK 

In this paper we explored the implications of the dis- 
persion relation of Eq. ^ on the ground state wavefunc- 
tion for several archetypical potentials. We found that 
the symmetry of the wavefunction could be changed by 
modifying the trapping potential. This feature would 
also be seen for more general dispersion relations where 
E{k) is non-quadratic and/or has multiple degenerate 
minima. Lin et al.'s recent realization of the disper- 
sion in Eq. ([2| offers an opportunity to experimentally 
test our predictions. 

The symmetry /antisymmetry of the ground state 
could be explored by either real-space probes (absorption 
imaging) or momentum-space probes (time-of- flight) . In 
particular, if l/fc^ is large compared to one's imaging 
resolution, one can simply count nodes or antinodes in 
the wavefunction. In Ref. 0, I/^l ~ 200nm, but this 
can be made longer by changing the angle between the 
Raman beams. 

An interesting use of time-of-flight would involve inter- 
fering outgoing waves with momentum k and — fc, giving 
a definitive measure of the symmetry/antisymmetry of 
the state. One should be able to map out the phase dia- 



gram seen in Fig.[6j changing /3 and ft to achieve different 
sized bumps in the physical potential and the dispersion 
relation. In particular, a clear transition should be seen 
going from a symmetric state localized in time and mo- 
mentum for f2 > AEl and /3 ^ 1, to an antisymmetric 
double- well state for /? ^ 1 and fl < 2El. This corre- 
sponds to a barrier whose depth is ^ 85 nK and a Raman 
coupling of order ~ 22 KHz. 

As illustrated in Fig. [ij for a symmetric potential, the 
transitions between different symmetry ground states are 
true crossings, and one cannot adiabatically change from 
one to another. However, if one introduces some asym- 
metry, these will become avoided crossings. Depending 
on details, adding interactions can either further smooth 
out these crossings, or sharpen them, leading to further 
hysteresis [Sj . Studying the role of interactions in these 
gases is an active area of research [7]. 
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